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Abstract 



The Clebsch method provides a unifying approach for deriving vari- 
ational principles for continuous and discrete dynamical systems where 
elements of a vector space are used to control dynamics on the cotangent 
bundle of a Lie group via a velocity map. This paper proves a reduc- 
tion theorem which states that the canonical variables on the Lie group 
can be eliminated, if and only if the velocity map is a Lie algebra action, 
thereby producing the Euler-Poincare (EP) equation for the vector space 
variables. In this case, the map from the canonical variables on the Lie 
group to the vector space is the standard momentum map defined us- 
ing the diamond operator. We apply the Clebsch method in examples of 
the rotating rigid body and the incompressible Euler equations. Along 
the way, we explain how singular solutions of the EP equation for the 
diffeomorphism group (EPDifT) arise as momentum maps in the Clebsch 
approach. In the case of finite dimensional Lie groups, the Clebsch vari- 
ational principle is discretised to produce a variational integrator for the 
dynamical system. We obtain a discrete map from which the variables 
on the cotangent bundle of a Lie group may be eliminated to produce a 
discrete EP equation for elements of the vector space. We give an inte- 
grator for the rotating rigid body as an example. We also briefly discuss 
how to discretise infinite-dimensional Clebsch systems, so as to produce 
conservative numerical methods for fluid dynamics. 

1 Introduction 

We are dealing with variational principles defined by an action (or cost, for the 
optimal control problem) 



whose Lagrangian (or cost functional) I : V t-^ W is defined on vectors ^ in a 
vector space V, subject to a condition imposed by a velocity map from the 
vector space V to the tangent space TqM of a manifold M at the point Q, 




(1) 



C^: V X M ^ TqM . 



(2) 
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The velocity map introduces the dynamics, 

Qit) = C^Qit), (3) 

where ^ G V and Q G TqM is tangent to the curve Q{t) in the manifold A/. 
Such variational principles arise in two different contexts: 

1. The optimal control context, in which one seeks solutions for Q{t) gov- 
erned by the dynamics Q that control the motion along a curve in an 
interval < t < T so as to minimise the cost in ^ for a given cost 
fimctional 

2. The Hamilton's principle context, in which stationarity SS = oi the 
action in ([T]) implies dynamical equations for ^ subject to the constraint 
imposed by the velocity map 

One approach that applies in both contexts was first introduced for ideal 
fluid dynamics in Serrin |Ser59| and in Seliger and Whitham [SW68j . This is the 
Clebsch approach for deriving variational principles for Eulerian fluid dynam- 
ics. A similar approach later emerged in the work of Bloch et al. |BCMR98] in 
the optimal control of rigid bodies. This approach enforces equation ^ through 
a Lagrange multiplier term in the action or cost. Doing so produces dynamical 
equations for Q and for the Lagrange multiplier P in terms of ^, together with 
a formula for ^ given in terms of Q and P. The Lagrange multiplier P is also 
the canonically conjugate momentum for Q in the corresponding Hamiltonian 
formulation, and the formula for ^ in terms of Q and P has special significance 
in the Hamiltonian framework. 

Such variational principles are said to be implicit Lagrangian systems 
and the subject has now reached a high state of mathematical development 
|YM06| . In this paper we take a "bare hands" approach to investigating this 
sort of problem. 

Section 2 describes the conditions under which the coordinate Q and its 
canonical momentum P in T* M may be eliminated from an implicit Lagrangian 
system in order to obtain a dynamical system for £ only. Answering this question 
summons a Lie algebra structure on V. That is, Q and P may be eliminated 
if and only if £^ corresponds to the action of some Lie algebra g on M; that 
is £| : g X M — s- TM for ^ G 0. In this case the dynamical system for ^ 
is always the Euler-Poincare equation for V with the appropriate ad* operator 
defined on the dual Lie algebra q* through the natural pairing induced by taking 
variational derivatives. The formula for ^ in terms of Q and P is then found 
via a cotangent-lift momentum map. This key result for the continuous case is 
stated and proved in Theorem [T^ 

In the case where V is the Lie algebra of vector fields X(M"), there is a 
choice of Lie algebra actions on different spaces e.g. Lie derivatives, left-action 
on embedding space etc. In each case the resulting system for £ is the Euler- 
Poincare equation for the diffeomorphism group. This is the EPDiff equation. 

Potential energy terms may also be introduced into the action in an implicit 
Lagrangian system by following standard procedure for Hamilton's principle. 
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This allows an extension of the Clebsch variational principle to obtain Euler- 
Poincare equations with advected quantities [HMR98| . Many equations of fluid 
dynamics may be obtained this way. Section 3 discusses the standard example of 
the incompressible Euler equations and provides references for other examples. 
Again all calculations are performed explicitly. 

Clebsch variational principles not only unify the subject, they also provide a 
systematic framework for deriving numerical integrators. A great deal of activity 
and rapid development of these variational integrators has recently transpired. 
See |BRM07| for an up-to-date survey of the subject, a bibliography and new 
results from the same viewpoint as the present paper. 

Section 4 discusses the potential for constructing numerical integration meth- 
ods by discretising the Clebsch variational approach in both space and time 
and deriving the resulting discrete equations of motion. The Clebsch approach 
provides a method of obtaining variational integrators |LMOW03] simply by 
discretising the Clebsch variational principle in both space and time. These 
integrators are symplectic and hence they fit into the backward-error analy- 
sis framework }LR05| . This means that they preserve the Hamiltonian within 
0{AtP) (p is the order of the method in time) and also preserve conservation 
laws associated with symmetries (provided that the symmetries are retained by 
the spatial discretisation). We show that for dynamics on finite-dimensional Lie 
groups, the discrete Clebsch variational principle results in equations where Q 
and P can once again be eliminated to produce a conservative integrator for the 
equation for ^ only. We give an example of a Clebsch integrator for the rigid 
body equation which after eliminating Q and P takes a particularly elegant 
form. The last section of the paper closes by discussing some possible direc- 
tions for applying the Clebsch approach and other recently developed parallel 
approaches for the case of infinite-dimensional systems. 

2 Clebsch variational principles 

Definition 1 (Velocity map). Suppose V is a vector space and M is a manifold. 
For each Q E M and ^ G we define the linear velocity map : V x M i-^ 
TqM. 

For a given element ^ of V, we use C to define velocity on TqM via 

Q = jC^Q, Q e M. 

Definition 2 (Clebsch action principle). For a given functional I : V -^M., the 
Clebsch action principle is 

s r iW)] + iP{t), Q{t) - ^mO) dt = 0, (4) 

Jtl ^ ' T* M 

where P is a Lagrange multiplier in TqM and {■,')t*m standard inner 

product on T* M . 
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Remark 3. The solutions of this action principle minimise J^^l{^)dt subject 
to the constraint that Q is directed by the velocity map C. The second term in 
^ imposes the definition of velocity by a Lagrange multiplier that will turn out 
to be the conjugate momentum in the Hamiltonian formulation. This will be 
shown in theorem ] 121 

To write down the general solutions, we need to define the diamond operator. 

Definition 4 (Diamond operator). Let C be a velocity map from V to M as 
defined above. The operator o : T* M — > V* satisfies 

{PoQ,^)y = -{PXiQ)r>M- 
where (■, ■)y is the inner product on V x V* . 

Remark 5. For the case where the velocity map is minus the Lie derivative, 
the diamond operation is the dual action of the Lie derivative. 

Remark 6. Later we shall see that for the case where the velocity map is a Lie 
algebra action, the quantity —P o Q is a cotangent-lift momentum map. 

Lemma 7 (Clebsch equations). 

The optimising solutions for the action principle ^ satisfy: 

51 „ ^ 

Q = c^Q, 

P ^ ^{TQC^fP, 



where one defines the variational derivative operations 
SI 



'^' = (77'^^) ' and {P,TqC^5Q) = {{TqC^YP,5Q) 



V 



Proof. 



= (5 [ ' I[i] + IP,Q- Cs^Q) dt, 

- r 

+ (^P,5Q - [TqC^) ■ 5Q - CsiO) ^^^dt, 

^(~P~{TQC^f P,5q) dt, 
and the result follows since (5^, 8P , and 5Q are arbitrary. □ □ 
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Lemma 8 (Legendre transform). Suppose I is chosen such that it is possible to 
solve for ^ from 51 /S^, i.e. there exists an operator G such that 

51 

Then the Clehsch equations for Q and P are canonical with Hamiltonian given 

by 

H={P,CG(PoQ)Q)-l[GiPoQ)]. 

Proof. The result follows directly from calculating the canonical equations for 
this Hamiltonian. □ □ 

Example 9 (The rigid body: Bloch et al. 1998). For example, consider the 
rigid body in coordinates {Q,Q) G TSO{3). The angular velocity of the rigid 
body is defined as 

Q^Qn with = - ne so(3). 

The Clebsch variational principle for the rigid body is given by 5S = with 
action S — J Ldt and implicit Lagrangian 

L{n, p, Q)^^n-m + p'^ -{Q- qn) 

for constant, symmetric I. The variations of L are given by: 

5 J Ldt ^ J 5n-{m- Q^P) + 5P^ ■ (Q - QVt) - 5Q^ ■ {P - PQ) 



dt. 



Thus, stationarity of this implicit variational principle implies a set of rigid 
body equations which first appeared in the work of Bloch et al. \BCMR98f on 
optimal control of rigid bodies 

m = -PoQ = Q^P, Q^CnQ = Qn, P = {TgCnQf P = Pn . 

In this particular example, one can study the dynamics of Q by calculating 

- Q^p + Q'^p, 

at 

= {QilfP + Q'^Pn, 

= n'^{Q^P) + iQ^P)Q, 

= fi^{m) + {m)n = [m, n]. 

Thus, the Clebsch equations yield the dynamics for Q by eliminating P and Q 
via the derivative of 51/55^. 

We will determine when it is possible to obtain a dynamical system for ^ by 
eliminating the Clebsch variables. First we need to make two more definitions. 
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Definition 10 (Closure). Let V he a vector space with C being a velocity map 
V X M ^ TM . The velocity map L is said to be closed if, for every pair of 
vectors u^v €z V , there exists a third vector w £ V such that 

CwQ = {{TqCv)Cu — {TqCu)C.u)Q, 

where TqCu is the tangent of Cu evaluated at Q. In that case, the velocity map 
induces an algebra structure on V with bracket 

which satisfies the Jacobi identity. 

Definition 11 (ad- and ad* -operators). Let C define an velocity map from V 
on M with induced bracket [•, •] : V x V ^ V . We define the ad-operator by 

adu V = —[u,v]. 

For m in V* we define the ad* -operator by 

(ad* m, v)y = {m, ad„ v)y . 

Theorem 12 (Elimination requires closure of the velocity map). 
Let V be a vector space, and let C be a velocity map from V x M ^ TM . Then 
the cotangent variables {P, Q) E T*M may be eliminated from the Clebsch 
equation if and only if the image of C in TM is closed under the Lie bracket. 
Furthermore, when this closure condition holds: 

1. C induces a Lie algebra structure on V , 

2. PoQ is a cotangent-lift momentum map, and 

3. ^ satisfies the Euler-Poincare equation: 

d SI SI 

d-tsi + '^^^rr'- 



Proof. First suppose that C is closed. For any vector w € V, 

[Definition of o] = -^(P,CwQ)t*m 
dt 

[P and Q equations] = {-{TqC^fP, C^Q)t'M + (-P, {TqC^)C^Q)t' m 
= {P,{-{TqC^)C^ + {TqC^)C^)Q)t'm 
[closure] = {P, -C[^^^]Q)t-m 
[Definition of o] = {P o Q, [|, w])v 

^{i),ad^w 
ad| 
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and hence 

d SI SI 

which is a closed system for ^. Hence, closure allows elimination of P and Q. 

Conversely, suppose that C is not closed. Then there exist u,v € V such 
that 

{TqCv)Cu — {Tq£u)£v 7^ £-w , (5) 

for any w G V. Now assume, aiming for a contradiction, that P and Q may be 
eliminated from the equations. For any Q, we may find a P such that 



61 , , 61 



= PoQ. 



Then 



= {P,{-{TqC^)£^ + {TqC^)£^)Q). (6) 

In order for this to be consistent, the left-hand side must be linear in £ and 
61/6^ (treating 61/6^ as a separate variable), and hence we may write 

d 61 r 

for some linear operator £u- Wc pair this with a vector, then use the Clcbsch 
variational equation and the definition of diamond to write 

d 51 \ /r , 

d75^^")'7v ^ \^"^(")'" 



Consequently, comparing with equation ^ yields 

^c^vQ ^ [(7q£u)^i; — {Tq£v)£u]Q 

for any Q, which contradicts the assumption ([5]). Hence, P and Q may be 
eliminated from the Clebsch equation if and only if the image of £ in TM is 
closed under the Lie bracket. 

The fact that — P o Q is a cotangent-lift momentum map comes straight 
from the definition: if £ is closed, then V is a Lie algebra with bracket induced 
by £, and the map £ is a Lie algebra action on M . The Hamiltonian for the 
corresponding cotangent-lifted action is 



= {P,£^Q)t' 



M, 
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and the momentum map J for this action is defined by the relation 

{J,i)v = {P,CuQ}t'M 
which matches Definition lU and hence 

J = PoQ. 

□ □ 

Remark 13. As a result of Theorem ] 1 'A solutions of the Clebsch equations may 
be composed by o into solutions of the following equation 



d SI „ 61 



whenever C is closed. This equation is called the Euler-Poincare equation and 
it describes geodesic motion on the space V whenever I is a quadratic functional 
(metric). 

To extend the scope of this framework it is useful to introduce potential 
energy terms to I in the action principle, i.e. let / be a function of the elements 
of M as well as V. 

Definition 14 (Clebsch action principle with potential energy). 
For a given functional I : V x AI ^ W, the Clebsch action principle is 

5 ' im). Q{t)] + {P{t). Qit) - C^it)Q)^^^/t = 0, (7) 

where P is a Lagrange multiplier in TqAI and (•, ■)x'm usual inner prod- 

uct on T*M . 

Lemma 15 (Clebsch equations with potential energy). The optimising solu- 
tions for the action principle are: 

SI „ ^ 

— = PoQ, 

Q = C^Q, 

P = + 



8 



Proof. 



P,SQ-{TQC^)-SQ-CsiQ) dt, 

/T'M 

^+PoQ,S^) +(5P,Q-C^q) 

04 /v ^ 'T'M 

-P-{TQC^f P+^,5q) dt, 

and the result follows since S^, SP, and 5Q are arbitrary. □ □ 

Theorem 16 (Elimination theorem with potential energy). Let V be a vector 
space, and let C be a velocity map from V — > TAl . The cotangent variables P 
and Q may be eliminated from the Clebsch equation with potential energy, if 
and only if the image of C in TM is closed under the Lie bracket. 

Furthermore, when the condition holds, ^ and Q satisfy the Euler-Poincare 
equation with advected quantities: 



^51 SI 51 ^ ^ 



Proof. The proof follows the proof of Theorem [121 using 
d 



dt 



-{TQCifP+^,C^Q) +{P,(TqCu,)C^Q 

□ □ 



T'M 



T'M 



3 Examples 

3.1 Singular solutions of EPDiff 

When V is the space of vector fields X(R"), under the conditions of Theorem 
1121 the Clebsch implicit variational principle yields the Eulcr-Poincare equation 
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for diffeomorphisms (EPDifF): 

d SI „ , SI . ,^ ,T 51 

One possible way to form a Clebsch principle for EPDiff is to consider the 
left-action of vector fields on the space of embeddings M = Emb(5, M") for 
some manifold S (such as the circle, or the sphere). For an embedding 

Q : 5 ^ M", 

the velocity map is defined by composition of functions 

£„Q(.s) = u{Q{s)), s e S. 
The Clebsch principle is then 

S f ' l[u] + {P,Q - u{Q))t'm dt ^ 0, 
Jti 

where the inner product in the second term is defined as 
^P(s,i). (Q(s,i)-M(Q(s,t)))ds. 

The diamond operator (o) is thus defined in this case by 

{{PoQ),u)v = {P^^QW'M 

= f P{s)-u{Q{s))ds, 
Js 

= / P{s) / 5{x~Q{s,t))u{x)dV{x)ds 

Js J M 

= J (^j P{s)5{x-Q{s,t))d^ ■u{x)dV{x). 

Consequently, one finds 

PoQ^ I P{s)S{x~Q{s,t))ds, 
Js 

which is the singular solution ansatz of [HM04| . 
One then calculates the Clebsch equations as 

A/ r 

— = PoQ= P{s,t)5{x-Q{s,t))ds, (8) 
Su Js 

Q = C^Q = u{Q), (9) 

P = -iTQC^fP^-J2Pk^ = -{Vu{Q)f-P, (10) 
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which we know to be canonicaUy Hamiltonian from Lemma [S] Furthermore, we 
know that Sl/Su satisfies the EPDifF equation from Theorem [T2l To verify this 
statement, take the inner product with a test function and differentiate in time, 
as follows. 

d / \ d ,„ _ , 



X{M) 



l^{P.w{Q)+P-^-Q)ds, 

[ - {{\7uf ■ P) ■ w{Q) + P ■ Vw{Q) ■ u{Q)ds, 
Js 

/ I - / Pd{x - Q{s))ds -Vuix) -wix), 
Jai \ Js 

+ J P5{x-Q{s))ds-yw{x)-u{x)^dVo\{x), 

^« / X(M) X*^" /X{M) 



This shows that 51 / 5u satisfies the weak form of EPDiff: 



dt\5u \ Su / j^^^i^ \5u / ^.(^^ 

Remark 17. In consonance with Theorem \lS\ and as discussed in \HM04^ the 
singular solution ansatz for EPDiff given by (0) is an equivariant momentum 
map arising from the cotangent lift of the action of vector fields corresponding 
to composition the left. 



3.2 Euler equations 

The Lagrangian for the incompressible Euler equations is 

f D 

l[u,D]= i — |Mp - i:))dVol(a;), 
Jm 2 

where D is the density, and p is a Lagrange multiplier introduced to enforce 
the constraint that the fluid is incompressible. Not unexpectedly, the quantity 
p turns out to be the pressure. 

There are several ways to write down a Clebsch variational principle for 
the incompressible Euler equations. However in order to obtain a set of vari- 
ables which include all possible solutions one needs to include at least d scalar 
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Lagrange multipliers where d is the spatial dimension of M. If, for example, 
we only use one Lagrange multiplier then only the vorticity-free solutions are 
obtained, as first noticed by Lin (see [CM87j for discussion and references). 

As described in |HK83| . one possible way to construct such a Clebsch vari- 
ational principle is to use the action of the diffcomorphisms on the Lagrangian 
labels, which satisfy 

ef = c^e"^ ^ - u-we"^ , £'^{x,o) = x'^, A = i,...,d. 

The Clebsch variational principle is then 

rt2 



5 f ^\u\' ~p{D~l)dYol{x) 
Jti Jm ^ 

+ I P ■ {£t+u-\/e)dVol{x) 

+ [ <p- {Dt + S/ ■ {uD))dVol{x)dt 
Jm 



= 0. 



Here we have also introduced an additional Lagrange multiplier (j) to enforce the 
dynamics of the density D which means we have d + 1 scalar Lagrange multipli- 
ers. The system could be reduced to d multipliers but this would unnecessarily 
complicate the exposition. 
The Clebsch equations arc: 

d£ 

f + V-(.I,) = 0, 
dP 

^ + V.(.P) = 0, 



SI 1 



Du^ -{S7ef -P- DVcj), 



dt ^ SD 2 

SI 



Su 



together with the constraint D = 1. After elimination one obtains 

^''af ^ + u ■ V{uD) + {Vuf ■ uD = ZJviiMp - Vj5, 
dD 

— + V-{uD) = 0, D^l 

which becomes 

du 

— + u- yu = -S7p, V • u = 0, 
ot 

after substituting D = 1. 

The Clebsch variational principle encompasses essentially all fluid theories. 
See |Hol02| for references and the example of a complex fluid. 
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4 Clebsch integrators 



In this section wc discuss the potential for constructing numerical integration 
methods by discretising the Clebsch variational principle in both space and time 
and thereby deriving the resulting discrete equations of motion. Any numer- 
ical method obtained this way will automatically be a variational integrator, 
and hence will inherit the accompanying conservative properties such as exact 
preservation of momenta where the discretisation preserves the symmetry of the 
continuous system, and the long-time approximate conservation of energy via 
backward-error analysis. See |LMOW03] for a review of variational integrator 
methods. 

In the finite dimensional case one simply needs to follow the variational 
integrator programme by finding a map which approximates the time-derivative, 
and then substituting it into the Clebsch variational principle. In the following 
sections we show how to do this for the case of finite dimensional Lie groups 
such as SO{Z). 

Definition 18. For a manifold M , define the derivative map (f>^t by 
(j)At : M X M -> TM, Q") ^ Tq„+iM, 

such that (/'At(Q"^\ Q") is an approximation to Q at Q" fLMOWOSf . 

Definition 19 (First-order discrete Clebsch action principle). 

For a given functional l{^), the first-order discrete Clebsch action principle is 

SF^sY, (liC) + (P", 0At(Q"+\ Q") - ieQ")) = 0. 

Lemma 20. The first-order discrete action principle given in Definition \19[ is 
optimised by P , Q and ^ satisfying the Discrete Clebsch equations 

^(0" + P"^Q" = 0, 

Proof. 



(<5P",0At(Q"+\Q")-ieQ") 



and the result follows after renumbering indices. □ □ 
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Remark 21. When M is a vector space, one may choose 

Qn+l _ fyi 



At 

and the equations become 

P"+i = P" - At (rQ„+i£^„+iQ)^ P"+\ 

Remark 22. These equations give the first-order symplectic method, known in 
jLROSf as symplectic Euler-A, for the Hamiltonian system given in Lemma\^ 
The adjoint method to this, known in JLR05^ as symplectic Euler-B, is obtained 
from the following discrete variational principle: 

N 
n=l 

Remark 23. Higher-order schemes may be obtained by replacing the first- 
order discretisation of the Q- equation enforced by the Lagrange multiplier P 
by Munthe-Kaas methods \MK98f (Runge-Kutta methods on Lie groups). 

4.1 First-order integrators for matrix groups 

In this section wc show how to construct a derivative map for the case of matrix 
groups. For this entire section: 

1. Q is a d-dimensional complex matrix. 

2. The velocity map acts by matrix multiplication by X from the right 

/:xQ = Qx. 

3. The inner product on TM is defined by the matrix trace using the transpose- 
conjugate operation (f) 

(P,Q) -Tr(PQt) =^Pfc,Qfc,, 

kl 

where the overbar indicates complex conjugation. 

4. The diamond operation is given by P o Q = — Q^^P since 

(PoQ,X) = -(P,£xQ> 

= -Tr(P(QX)t) 
= - Tr (Qtpxl') 

= -(Q^P,x). 
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We first define the exponential map when the velocity map £x is defined by 
right matrix multiplication. 

Definition 24 (exponential map). The exponential map corresponding to the 
velocity map Cx is the solution of the equation 

^ exp(tX) = Cx exp(<X) := exp(iX)X, exp(O) = Id . 

We also define the logarithm map: 

Definition 25 (logarithm map). Let C be a velocity map. The logarithm map 
is defined by 

log(exp(X)) = X, 

for any vector element X of V . 

For a pair of group elements Qi, Q2; we seek a vector X such that 

Q2 = Qiexp(A«). 

This allows one to use £xQ2 as an approximation for the time derivative at Q2, 
thereby motivating the following definition. 

Definition 26 (Discrete approximation of time derivative). 

For two group elements Qi, Q2 an approximation to the time derivative at Qi 
of a solution which passes between Qi and Q2 in time At is 

</'At(Q2,Qi) = Qi^l^(Qr'Q2), 

where log : V is given by 

l^iA) = log(A) + OiiA ~ If), 

for some positive integer p. 

Let us now construct the Clebsch integrator. 

Definition 27 (First-order discretisation). We replace the time integral in the 
Clebsch variational principle by a sum and substitute our discrete approximation 
of the time derivative to get 

^Zj ('(X""') + {^p""\^Q""'i^((Q""')"'q") -Q"~'x""'^) =0. 

Note that the inner product is taken at the point Q"~^. 
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Theorem 28. The Clebsch variational principle given in Definition \27\ leads to 
the following first-order symplectic Euler discretisation: 

^(X"-i) = (Q"-i)tp"-i, (11) 
Q" = Q"-iexp(AtX"-i), (12) 
P" = ((Q")y'((r53^(AtX")l^)y'(exp(AtX"-i))^ 

(ikp(Atx-i)l^^ (Q"-i)^ P"-i (exp(AiX"))-^ ,(13) 
where cxp is the inverse of the log operation. 

Remark 29. The discrete Clebsch momentum map takes the expected form. 
Proof. The variational principle becomes 

TV 

,(5X^ 



^(5Q" -(5Q"-i (Q"-i)"^Q" 



The discrete Clebsch equations are then 

= |^(X"-i) - (Q"-^)tp"-\ (14) 

OA 



= Q"-ilog((Q"-i) 'Q") - AtQ"-iX"-\ (15) 

= 



P"((Q")"'Q"+i) + (16) 
^bi((Q")^^Q"+i)Q-X"yP". (17) 



After rearrangement, equation (|15p becomes 

Q"+i = Q"cxp(A<X"), 

16 



and the last equation simplifies to 

= ((Q"-i)-^Q")'(T(Q„_,)-.Q„bi)'(Q"-i)^P"-i + 

(r(Q„)-iQ„,,bi)VQ")tp" ((Q")-^Q"+i) , 

= (cxp(A«"-i))^ (lkp(A«"-Ml^)^ (Q""')^ P""' - 

fe(A«")bi)^ (Q")^ P" (cxp(AtX")) , 

as required in the statement of Theorem [28l □ □ 

Corollary 30. P and Q can be eliminated from the equations arising from the 
discrete variational principle in Definition \27\ to obtain the equation 

= (cxp(AiX"-i))^(r^(A«"-i)l^^^(X"-i)- 

(7kp(AtX")l^)^ ^(X") (cxp(AtX")) . 

Proof. Substitute equation (fTTj) into equation (fT2|) . □ □ 

Cayley transform methods In the following example, we derive an integra- 
tor for the rigid body equations by approximating the exponential map using the 
Cayley transform, which preserves the property of mapping from the Lie algebra 
into the group. This property of the Cayley transform has long been used for 
ensuring that numerical schemes preserve Lie group structure. In [AK W93] . the 
Cayley transform was used to reconstruct the attitude of a rotating rigid body 
from numerical solutions of the body angular momentum equation obtained 
from the midpoint rule. It was noted that the conservation of the Casimir 
||m|p (where m is the angular momentum) was necessary to obtain conserva- 
tion of spatial angular momentum. [LS94j proposed to transform a Hamiltonian 
system on a Lie group onto the Lie algebra using the Cayley transform (rather 
than the exponential map), and integrating the resulting equation numerically. 
This approach was developed in [IseOlj which suggested that, rather than in- 
tegrating the Lie algebra equation using a Runge-Kutta method (producing a 
Cayley Munthe-Kaas method [MK98| ). one could obtain an efhcient scheme by 
using a truncated Magnus expansion and, if a suitable numerical quadrature is 
used to approximate the integrals in the series, one obtains a time-reversible 
method of even order. In the example below, we embed the Cayley transform 
into the discrete Clebsch variational principle; higher-order schemes could be 
produced by using the methods of Munthe-Kaas and Iserles. 

Example 31 (Rigid body integrator). For the case where V is so (3) and acts 
on 50(3), we may approximate the exponential map to first-order using the 
Cayley transform 
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We obtain a corresponding approximation to the logarithm by taking the inverse: 

logA = cxp^^A, 

= 2(A - + = (A -/)(^/+^) \ 

= {A-I) + 0{{A-I)^), 
= T^{A) + 0{{A-ir). 

Our approximation to the time derivative is then 

0At(Q"+\Q") = Q"bi((Q")-iQ"+i), 

which is the usual linear difference with a projection applied to ensure that (/>At 
is in 5o(3). 

The Q- and H.- equations are then 

^(X") = (Q")^P", (18) 

/ X" \ / X" \ 

Q«+W/_At— j = QM/ + At— j. (19) 



The (^-component of the variational principle, which gives rise to the ¥ -equation, 
is 

N 
n=l 

Making use of the formula 

toi((Q"-^)-^Q") '^^^Q"'P''^"^ = 5((Q"-^)-^Q"), 
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we 



have 



^(log((Q"-i)-iQ")) = (/-iog((Q"-i)-iQ" 



2 



'5((Q"-')"'Q") 



-1 



'((Q"'^)'^Q")+/^ 



2 

-1 

-vn-l\-l 



(Q^ 



2 

((5Q"-<5Q"-i(Q"-i)-iQ 
Consequently the V -equation in this formulation is 



2 

((Q")-iQ"+i) +/' 



■((Q")-'Q 



n+l\ 



(Q")^P"((Q")-^Q"+i)^ 

Making use of equations hl^l9fl allows this to be rearranged into the more com- 
pact form, 



(exp(A<X"-i)) 

(5X I \ 2 

oX \ V 2 
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Finally, making use of the Cayley transform approximation, 

exp(AiX) =(^I + Af|^ (^I - At| 
and its consequenee, 

X 



/ + exp(AiX) = 2 / - At 



we obtain the reduced equation for X; 
This finally rearranges to become 

|(x") ^ ^(x"-V^(^(x"-)x"- + |(x")x^ 

+ ^ ((X"-i)^^(X"-i) + (X"f ^(X") 

(AO' ^ x"i^(X")X" -x"-ii^(x"-i)x"-i 

which is the discrete rigid body equation obtained from this choice of discrete 
Clebsch variational principle. 

Remark 32. The Clebsch integrator for the rigid body obtained using the Cayley 
transform approximation for exp is equivalent to the CAY-integrator obtained 
from the discrete Hamilton-Pontryagin principle in JBRMO^. In that case, the 
equations are obtained by extremising a functional on a Lie algebra (in this case 
the kinetic energy as a function of the body angular velocity X ) subject to the 
constraint that QQ~^ = X; in the case of the CAY-integrator the exponential 
map is again discretised using the Cayley transform. The Hamilton-Pontryagin 
principle provides an alternative viewpoint to the Clebsch principle; the extra 
feature in the Clebsch framework is the role of the o- operator as a momentum 
map. 

A plot of the dynamics obtained from this discrete integrator is given in 
Figure [TJ 



5 Summary and outlook 

This paper has discussed Clebsch variational principles from the point of view 
of a velocity map Q = C^Q which allows the dynamics Q{t) on a manifold to 
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Figure 1: A plot showing numerical integration over 100 periodic orbits using the 
discrete Clebsch integrator for the rotating rigid body with moment of inertia 
eigenvalues (0.5,0.6,1) and At = 0.1. The good energy conservation and exact 
angular momentum conservation are illustrated through the persistence of the 
periodic orbit structure over this long time integration interval. 

be controlled by a time-series ^{t) of elements of a vector space, using Lagrange 
multipliers P{t). Theorem [T^ shows that Q and P may be eliminated from the 
resulting Clebsch equations, if and only if the velocity map is a Lie algebra 
action on Q £ M. The Clebsch framework for velocity maps thus has a clear 
connection with the theory of Euler-Poincare reduction; namely, the equations 
obtained are the Euler-Poincare equations on the dual of the Lie algebra. For the 
continuous time case where the velocity map is assumed to be a Lie derivative, 
this connection is not unexpected. 

Examples in the paper included the finite-dimensional rigid body equation, 
and two infinite-dimensional examples: the singular solutions of the EPDiff 
equation and the incompressible Euler equation. In the EPDiff example, the 
Clebsch framework derives the singular solutions as a family of momentum maps. 

Finally the paper showed how discrete Clebsch variational principles can 
be used to produce numerical methods for Clebsch equations. For the case 
of finite-dimensional Lie groups, in which the variational principle need only 
be discretised in time, one may again eliminate Q and P using the discrete 
approximation for the time derivative in Definition [55] to obtain a conservative 
numerical method in terms of ^ only. The example of discretisation for the 
rigid body resulted in the CAY integrator for the associated EP equation in 
|BRM07j . Possible extensions of this technique would be to obtain higher- 
order time-integration methods based on Rungc-Kutta/Munthc-Kaas methods 
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|MK98| or Magnus methods [IseOlj . and to apply the Clebsch mtegrator to 
systems with potentials such as the heavy top. 

In the case of infinite-dimensional systems, it is necessary to discretise the 
variational principle in space as well as time, if one wishes to eliminate Q and 
P thereby obtaining a closed discrete equation for ^, Theorem requires the 
spatial discretisation of the velocity map to remain a Lie algebra action on the 
discretised space. This is a key step for making future progress in applying the 
Clebsch method for discretising fluid dynamics and other infinite-dimensional 
evolutionary systems. In the special case where the Lagrangian is at most linear 
in space-time derivatives (without higher derivatives) , then the resulting PDE is 
multisymplectic (see jBROlj and cited papers). Any discrete Clebsch variational 
principle leads to a multisymplectic Clebsch integrator. See |CHH07j for more 
details. 
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